Nonsequential Two-Photon Double Ionization of Atoms: Identifying the Mechanism 
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We develop an approximate model for the process of direct (nonsequential) two-photon double ionization of 
atoms. Employing the model, we calculate (generalized) total cross sections as well as energy-resolved differ- 
ential cross sections of helium for photon energies ranging from 39 to 54 eV. A comparison with results of ab 
initio calculations reveals that the agreement is at a quantitative level. We thus demonstrate that this complex 
ionization process is fully described by the simple model, providing insight into the underlying physical mech- 
anism. Finally, we use the model to calculate generalized cross sections for the two-photon double ionization of 
neon in the nonsequential regime. 
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PACS numbers: 32.80.Rm, 32.80.Fb, 42.50.Hz 

Correlated dynamical processes in nature poses unique 
challenges to experiments and theory. A prime example of 
this is the double ionization of helium by one-photon im- 
pact, which has been studied for more than 40 years. How- 
ever, it is only during the last 15 years or so, that advances 
in theory, modeling and experiment have enabled scientists 
to gain a deeper insight into the role of electron correlations 
in this ionization process The corresponding prob- 

lem of two-photon double ionization of helium, in the photon 
energy interval between 39.4 and 54.4 eV, is an outstanding 
quantum mechanical problem that has been, and still is, sub- 
ject to intense research worldwide, both theoretically ifs Uigll 
and experimentally, employing state-of-the-art high-order har- 
monic and free-electron (FEL) light sources [i23Ll24ll . 
Despite all the interest and efforts that have been put into this 
research, major fundamental issues remain unresolved. What 
characterizes this particular three-body breakup process is that 
the electron correlation is a prerequisite for the process to oc- 
cur, i.e., it depends upon the exchange of energy between the 
outgoing electrons, and as such it represents a clear departure 
from an independent-particle picture. 

In this Letter, we present a novel approximate model for the 
direct or nonsequential two-photon double ionization process 
in helium, sketched in Fig. [T] (a). We show that the simple 
model predicts the essential features of the process, even at 
a quantitative level, which is quite surprising given the very 
high complexity of the problem. In particular, we find very 
good agreement between the model predictions and the results 
obtained by solving the time-dependent Schrodinger equation 
from first principles, regarding (generalized) total cross sec- 
tions as well as energy-resolved differential cross sections for 
the process. The proposed model may be generalized to ac- 
count for direct double ionization processes in multi-electron 
atoms. We demonstrate this by calculating the generalized 
cross section for nonsequential two-photon double ionization 
of neon. 

Few-photon multiple ionization of noble gases beyond he- 
lium have been studied experimentally in some detail i23l - l26ll . 
but to the best of our knowledge, the cross section for the 
nonsequential two-photon double ionization process has not 
yet been obtained. Therefore, we hope that our results will 



encourage further investigation of nonsequential double ion- 
ization processes in various noble gases. 

Reducing a complex quantum mechanical problem to a 
simple and transparent model problem, while retaining the 
essential physics, is very useful in order to access the un- 
derlying physics li, 27, 28]. With such a goal in mind, we 
will now outline a possible physical mechanism for the non- 
sequential two-photon double ionization process in an atom, 
and then proceed to construct a simple quantum mechanical 
model which implements these ideas. The idea behind the 
model is that the electrons are considered to be distinguish- 
able particles that can absorb one photon each. However, in 
order to include the effect of the first emitted electron on the 
second one, we impose the additional but important constraint 
that the absorption of the second photon, by the second elec- 
tron, can only occur after the first photon absorption. In this 
way, and according to the principle of conservation of energy, 
the first electron may transfer energy to the second electron as 
it is emitted, allowing for the nonsequential ionization process 
to take place. 

The starting point of our model is the single-active electron 
approximation (SAE) where both electrons are considered to 
be independent particles and treated differently in that they are 
both assumed to move in their respective ionization potentials. 
That is, the 'outer' electron moves in an effective potential set 
up by the nucleus of charge Ze {e is the elementary charge), 
the 'inner' electron and the Z — 2 other electrons. The inner 
electron sees a corresponding screened potential given by the 
nucleus and the Z — 2 remaining electrons. We will label these 
two different cases simply by 'A' and 'B', respectively. Fol- 
lowing this procedure, the wave function of the ground state 
may be approximated by the product ansatz 



'P(rA,rB) = \ifA{rA)WB{rB) 



(1) 



where \//4 and y/s refer to the one-electron wave function of 
electron A and B, respectively. 

Now, the first ionization event in the direct two-photon dou- 
ble ionization process can be represented by the one-electron 
dipole coupling between the ground state wave function of ei- 
ther A or B, i.e., the state jA,/?^) or Iz?,/?]]), and their respec- 
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FIG. 1. (color online), a) Sketch of the direct two-photon double 
ionization process in helium. The abbreviation SI and DI stands for 
single and double ionization continuum, respectively, whereas the ar- 
rows illustrate the photons that are absorbed by the system, b) Sketch 
of the model process for two-photon double ionization (see text for 
details), c) Matrix representation of the model Hamiltonian, for the 
case where the outer electron is emitted before the inner electron (see 
text for more details). Atomic units (a.u.) are used in the figure (1 
a.u. of energy corresponds to 27.2 eV). 



tive continuum states, \A,Ea) and \B,Eb}, where and Eg, 
and Ea and Eb represent the energies of the ground and con- 
tinuum states, respectively. In the product basis representa- 
tion ([TJ, with the length gauge formulation of the light-matter 
interaction, the dipole coupling matrix elements may be writ- 
ten on the following simple form (for the case where electron 
A is assumed to be emitted first). 



{A,E^\~eEit) ■rA\A,EA) {B,E^s\B,E^s) . 



(2) 



where E(f ) is the time-dependent electric field that defines the 
laser pulse, which is assumed to be Unearly polarized along 
the z-axis. These coupling elements are related to the one- 
photon (one-electron) photoionization cross section via the re- 
lation S 

cja = 4K^a (£a I {A,E^\za\A,Ea)\^ , (3) 
where a is the fine structure constant. 
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FIG. 2. (color online). Total integrated (generalized) cross section 
for the nonsequential two-photon double ionization of helium. Black 
line: present model result obtained with Eq. l[5](; open (blue) circles: 
ab initio result of Feist etal. |T 5] obtained with a 4 f s pulse ; and open 
(red) squares: corresponding ab initio result of Nepstad et al. 11^ . 
The vertical lines define the two-photon direct double ionization re- 
gion. 



At the instant of ionization of electron A, electron B remains 
unaffected. However, once electron A has absorbed its photon, 
we allow for the possibility that electron B (but not A) can be 
hit by a second photon. This secondary process is included 
into the model by introducing additional dipole couplings be- 
tween the B ground state and its corresponding one-electron 
continuum states in the following way: 



{B,E°s\-eE{t)-rB\B,EB)8iEA,E'A) 



(4) 



Note here that there are only non-vanishing couplings be- 
tween SAE states (system A) of the same energy, i.e., the 
resulting coupling matrix attains a very simple structure, as 
shown in Fig. [T] with typically only a few hundred different 
couplings. The same procedure may also be followed with A 
and B interchanged, however, this will neccessarily yield the 
same result, and therefore need not explicitly be considered. 

The couplings (O and ^ and the mentioned constraints, 
along with the corresponding diagonal energies, constitute the 
entire model that we propose. To this end, we would like to 
add that all excited, bound states have been left out of the 
model, as they play no role in the present context. As a matter 
of fact, despite the extremely simple form of the model matrix 
elements, with no explicit presence of the coiTelation poten- 
tial, it actually allows for the possibility that the two electrons 
exchange energy in the excitation process. Thus, both elec- 
trons may be emitted into the continuum even though the en- 
ergy of the secondary photon may not itself be sufficient to 
eject the inner electron into the continuum. 

Applying second order perturbation theory to the resulting 
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FIG. 3. (color online). Electron energy distribution for two-photon 
double ionization of helium at photon energies of 44.9 and 51.7 eV. 
Solid (black) line: model result; and dashed (red) line: ab initio re- 
sult. 



model Hamiltonian, one can show that the single-differential 
cross section for the direct two-photon double ionization of an 
atom is simply given by 



dE 2^-'^ ' 



-f{2fico+E^ 



(5) 
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where Oa and Ob now refer to the total one-photon single ion- 
ization cross sections of A and B, respectively, E is the ex- 
cess energy, and where we have explicitly accounted for the 
exchange symmetry of identical particles and the possibility 
that either the inner or the outer electron is emitted first. At 
this point we would like to emphasize that the only parame- 
ters needed in order to calculate the nonsequential two-photon 
double ionization cross section within the model framework, 
is the effective binding energies of electron A and B, as well 
as their respective one-photon single ionization cross sections. 
For instance, for helium all these parameters are well known. 
The model may straightforwardly be generalized to account 
for e.g. nonsequential three-photon triple ionization processes 
in atoms. A more detailed exposition of the model and a 
derivation of the perturbation theory expression for the cross 
section, will be outlined in a forthcoming communication. 

In Fig. |2] we compare the total cross section obtained us- 
ing the approximate model, Eq. (|5]l, (black line in the figure), 
with the ab initio result of Feist et al. BlSll (blue circles) and 
Nepstad et al. 1 19] (red squares), both of which were obtained 
by solving the time-dependent Schrodinger equation of he- 
lium from first principles. The model result is obtained using 
tabulated values for the absolute one-photon photoionization 
cross section of helium, as obtained experimentally by West 



and Marr fso']. Figure[3]shows corresponding energy -resolved 
single-differential cross sections at two selected photon ener- 
gies, 44.9 and 51.7 eV. As a matter of fact, the agreement 
between the model result and the ab initio results is almost 
perfect in Figs. |2] and [3] in particular for the lower photon 
energies, demonstrating the strength of this extremely simple 
model in predicting accurate values for the generalized cross 
section in direct two-photon double ionization processes. For- 
mula (|5]l predicts a sharp rise of the total cross section in the 
vicinity of the threshold at 54.4 eV, which is in agreement with 
recent ab initio calculations lfH[l3i[l5i [l8i[l9il3lll. 

As mentioned in the introduction, the problem of nonse- 
quential two-photon double ionization of helium has been 
subject of intense research in recent years, and accurate pre- 
dictions for the generalized cross section remain elusive fsj- 
|24|1 as the values obtained for the cross section for the reac- 
tion may differ by as much as an order of magnitude. On the 
theoretical side, the great discrepancies that remain between 
different approaches are usually ascribed to the different ways 
electron correlations are handled in the final state. To this end, 
we hope that the predictions of the present model study may 
shed new light on this controversy. 

Having justified the validity of our simple approach, we 
now turn to a more complex problem, namely the process of 
nonsequential two-photon double ionization of neon. Insert- 
ing, in Eq. (|5]l, the coiTect first and second ionization ener- 
gies of neon, i.e., 21.6 and 40.9 eV, as well as experimental 
values for the photoionization cross sections of Ne [30] and 
Ne+ llsill . obtained using synchrotron radiation, the result- 
ing model prediction for the double ionization cross section 
is shown in Fig. |4] (upper panel). The lower panel shows the 
corresponding electron energy distribution at three selected 
photon energies. Interestingly, at lower photon energies, the 
energy distribution exhibits a maximum (negative concavity) 
when both electrons are emitted with the same energy, while at 
higher photon energies the distribution is U-shaped. In sharp 
contrast to this trend, for helium, the model yields a U-shaped 
energy distribution for all photon energies (see Fig. [3]). 

In conclusion, we have implemented an approximate and 
very simple model to study the two-photon double ionization 
process of helium in the direct regime, i.e., at photon ener- 
gies below 54.4 eV where the sequential ionization process is 
energetically inaccessible. We have investigated the validity 
of the model by calculating generalized total cross sections 
and energy-resolved differential cross sections and compared 
the model results with corresponding results obtained by ac- 
curate ab initio calculations. Quantitative agreement between 
model results and the full results was achieved in all consid- 
ered cases, demonstrating the general validity of the model for 
the two-photon double ionization process. Finally, we have 
obtained the cross section for nonsequential two-photon dou- 
ble ionization of neon, demonstrating that the model has a 
great potential to be used in studies of nonsequential multi- 
photon multiple ionization processes in more complex atomic 
systems. This is an avenue of research we plan to pursue in 
the future. 
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FIG. 4. Upper panel: generalized total cross section for the process 
of two-photon double ionization of neon in the direct regime. The 
model result is obtained using Eq. (|5}, inserting available experimen- 
tal values for the total one-photon single ionization cross sections of 
neon |30] and Ne+ |32], respectively. The vertical lines define the 
two-photon direct double ionization region. Lower panel: normal- 
ized energy distributions at various photon energies. 
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